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ABSTRACT 

The interstellar scattering responsible for pulsar parabolic arcs, and for intra-day vari- 
ability of compact radio quasars, is highly anisotropic in some cases. We numerically 
simulate these observed phenomena using totally anisotropic, power-law models for 
the electron density fluctuations which cause the scattering. By comparing our results 
to the scattered image of PSR B0834-t-06 and, independently, to dual-frequency light 
curves of the quasar PKS1257-326, we constrain the nature of the scattering media on 
these lines of sight. We find that models with spectral indices slightly below /3 — 3, 
including the one-dimensional Kolmogorov model, are broadly consistent with both 
data sets. We confirm that a single physical model suffices for both sources, with the 
scattering medium simply being more distant in the case of B0834-I-06. This reinforces 
the idea that intra-day variability and parabolic arcs have a common cause in a type 
of interstellar structure which, though obscure, is commonplace. However, the implied 
gas pressure fluctuations are large compared to typical interstellar pressures, and the 
magnetic stresses are much larger still. Thus while these scattering media may be 
commonplace, their underlying dynamics appear quite extraordinary. 

Key words: Scattering - turbulence - ISM: general - ISM: structure - pulsars: 
individual: B0834-h06 - BL Lacertae objects: individual: PKS 1257-326 



1 INTRODUCTION 

Radio waves are scattered by inhomogeneities in the ionised 
interstellar medium (ISM). The appearance of astronomical 
objects is thus altered and the ISM leaves its imprint in 
the received radiation. The effects of propagation include 
angular broadening and wander of images, variation in the 
arrival time of impulses, and modulations of the flux density 
of compact radio sources (i.e. scintillations). These various 
effects provide powerful tools with which to study the ISM, 
including exploration of small-scale features that cannot be 
seen with any other methods (see, for example, the various 
contributions to|Haverkorn fc Goss|2007 l. 



Although many aspects of the available data can be 
successfully described using a model based on distributed 
Kolmogorov turbulence in the ISM, there are also some sig- 



nificant departures from this picture (Rickett 19901. Per- 



haps the most striking examples are the Extreme Scatter- 
ing Events ( [Fiedler et al]|1987[ [l994| , which appear to re- 
quire discrete, high-density plasma concentrations of size 
~ lO^** cm. And there are two other phenomena which 
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are indicative of highly localised scattering regions dom- 
inating the rest of the line-of-sight, namely the large- 
amplitude intra-day variations (IDV) of some compact radio 
quasars (Ked ziora-Chudczer et al.|1997| [Dennett-Thorpe &:| 
|de Bruyn||2000 Bignall et al.[[2003 l, and the parabolic arcs 
seen in the power-spectra of the dynamic spectra of radio 
pulsars | [Stinebring et al.|2001[ ). At present it is unclear how 
these three phenomena relate to the broader ISM, or how 
much they contribute to the scattering on a "typical" line- 
of-sight. 

There is strong evidence for high levels of anisotropy in 
the scattering material responsible for IDV. These indica- 
tions come from the quasi-sinusoidal character of the light- 
curves, and from modelling of the annual variation in scintil- 
lation time-scales, along with the scintillation pattern arrival 
time delays between widely separated telescopes (Rickett, 
Kedziora-Chudczer & Jauncey 2002; [Dennett-Thorpe fc de| 
[Bruyn[|2003| [BignaU et al.[[2006[ ). The anisotropy, as mea- 
sured by the major- to minor-axis ratio. A, must he A ^ 1, 
and it is possible that A ^ 10. In this circumstance to- 
tally anisotropic, one-dimensional scattering models are of 
interest as a simple, limiting case which can provide a good 
approximation to the actual circumstance. It was demon- 
strated by Walker, de Bruyn and Bignall (2009) that a one- 
dimensional {A — )■ oo) kinematic model offers a good ex- 
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planation for the two-station time-delays and annual cycle 
of the scintillation time-scale of the very rapid IDV quasar 
J1819-h3845. 

In the case of pulsar parabolic arcs, highly anisotropic 
scattering models are known to give arcs which are very 
sharply defined ( jCordes et al.|2006|| Walker et al.|2004l ). And 
a quasi-one-dimensional interpretation of the pulsar data is 
supported by observational studies which are directly sensi- 



media is far from settled, either observationally (Blandford 



five to the angular structure of the scattered image ( Walker 
et al.|2008||Brisken et al.|2010[ ) 



In the present paper we explore the constraints that can 
be placed on the parameters of a one-dimensional, power law 
model of electron density fluctuations, in the context of IDV 
and parabolic arcs. We do so using extensive numerical sim- 
ulations, which are compared against observations — visu- 
ally, and using summary statistics. We prefer this approach 
to any comparison based on ensemble-average models be- 
cause the data are acquired over short intervals, which are 



far from the ensemble-average regime ( Goodman & Narayan 



1989[ ). 

Section [2] introduces the model and simulation tech- 
nique. Section [3] derives model constraints from the observed 
image of a highly scattered pulsar, B0834+06, and sec- 
tion |4] formulates model constraints from the dual-frequency 
light curves of the IDV blazar PKS1257-326. Section 5 sum- 
marises our constraints, and in section 6 we draw inferences 
concerning the nature of the scattering media. 



2 SCATTERING MODEL AND SIMULATIONS 



We assume that there is a single thin screen (Bramley 19541, 
at an effective distance Dcff, that dominates the scattering, 
as appropriate for both objects of interest ( [Hill et aL]|2005[ 
Bignall et al.|2006 l. We model the electron column density, 
Ke, fluctuations as Gaussian random variables which are spa- 
tially correlated, with the power spectrum of Ke being ( |Lee 
fc Jokiph|19"76| [Narayan fc Goodman] 1989[ ) 

-Q/Qi 
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The spectrum of phase fluctuation (f> = Xr^He follows triv- 
ially. In equation , A is the wavelength of observations, 
is the classical electron radius and If ~ \Dcs/2tt is the 
Fresnel length. The parameters a, /?, Qo and Qi describe the 
spectrum itself. The quantity a is, up to a factor 0(1), the 
square of the characteristic phase variation at a transverse 
separation equal to the Fresnel scale and thus it gauges the 
amplitude of scintillations (e.g. |Narayaii||1992| ): a ^ 1 in- 
dicates weak scintillations, whereas a S> 1 corresponds to 
order-unity flux variations. Consequently we refer to a as 
the scattering strength. At sufficiently strong scattering, for 
power- law spectra with vanishing inner scale, a ~ it^^^ in 
terms of the ur parameter defined by equation 2.5 of |Rickett| 
IT990I ). 

Spectra with slopes /3 < 3 will be referred to as "shal- 
low" while those with /3 > 3 as "steep" . The classical theory 
of turbulent cascades ( |Kolmogorov| [T94T| predicts a shal- 
low slope oi P = 8/3, in the one-dimensional case (i.e. two- 
dimensional, but projected onto the plane of the sky), for the 
inertial range of a turbulent cascade. We emphasise, how- 
ever, that the physics appropriate to interstellar scattering 
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1 or theoretically (Iroshnikov"1963 
h & Sridhar|1995, Boldyrev,2006y 



And this point is underlined by our current lack of under- 
standing of the particular phenomena under consideration 
here. In this paper we consider a range of both shallow and 
steep spectra, 2 < /3 < 5. 

The outer scale Qq^ corresponds roughly to the largest 
scales at which the phase fluctuations remain correlated, and 
this scale dominates the statistics of phase, especially for 
steep slopes. Estimates for the value of Qo on different lines- 
of-sight vary consi derably (|Rickett|1990| Ar mstrong, Rickett 
& Spangler 1995, Chepurnov & Lazarian 20101. Since the 



outer scale cannot exceed the size of individual strongly- 
scattering plasma patches, intermittency of the IDV phe- 
nomenon sets an upper limit on the outer scale of the me- 
dia responsible for that phenomenon. Intermittency on a 
time-scale of months was seen in the intra-hour variable 



PKS0405-385 ( Kedziora-Chudczer 20061. Similarly, after a 



decade of rapid variations J1819-f 3845 ceased to scintillate 
( Cimo|2008| ). Finally, amongst a sample of hundreds of com- 
pact radio sources when IDV was seen it was often intermit- 
tent on a timescale of months (Lovell et al. 20081. These 



observations indicate that Q'^^ < 10^* cm. 

The inner scale, , is that of the energy dissipation 
in the classical theory. The microscopic processes which may 
be responsible for such dissipation include ion-neutral colli- 
sions and Landau damping ( Spangler|1991 1 . The inner scale 
affects strongly the phase derivatives and small-scale cor- 
relation properties of the screen ( jSpangler fc Gwinn|p90| ). 
A recent estimate, based on the impulse-response function 
for strongly-scattered pulsars ( Rickett et al.|[2009 l puts the 
inner scale below 10^ cm alt hough similar co nsiderations 
yielded values close to 10* cm ( Bhat et al.|2004" l earlier, and 
inner scales as large as 10^^ cm have previously been sug- 
gested ( Coles et al.|1987 l. To date there has been no deter- 
mination of the inner scale specific to the highly anisotropic 
scattering media responsible for IDV or parabolic arcs. 

We simulate realisations of the electric field envelope. 



w, in the manner of Martin & Flatte (1988); Narayan & 



Goodman] ([19891 . In modelling PSR B0834-I-06, a large spa- 
tial dynamic range A'^ > 10^ is needed in order to keep 
edge effects from corrupting the simulated observations. For 
one-dimensional simulations, this is feasible, which allowed 
us to repeat the simulations on demand, build statistically 
meaningful samples, and explore the parameter space in an 
adaptive manner. When extracting light curves, we com- 
puted the squared magnitude of the field at the sampling 
point, F„ = \uf. For IDV sources, incoherent averaging 
was performed as well - i.e. smoothing out the flux den- 
sities Fu{R) - with a Gaussian source profile. When statis- 
tical measures of the light curves were used for comparing 
data to the simulations - e.g. auto- and cross-correlation 
functions (ACF and CCF) - procedures were applied to the 
simulated flux densities which exactly paralleled the treat- 
ment of the comparison data. The angular structure of the 
images was estimated by Fourier transforming the spatial 
sampling of the field, so as to match the way the single dish 
data were analysed, then calculating its squared magnitude, 
1(6) = \j Axu^{x)e:^p(2-nivxe/c)\^. 

Altogether, we have five parameters in the model: 
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Figure 1. The scattered image intensity 7 as a function of angular 
distance 8 of PSR B0834+06, determined from the delay-Doppler 
image of Walker et al (2008). The low-intensity baseline of this 
image, at / ~ 10^*, is just the instrumental noise-floor. 



strength a, slope /3, outer and inner scales Qo and Qi and 
the source intrinsic size ae - or, equivalently, its brightness 
temperature Tb, as per = 2-jTagkTb/X'^ ■ Exploring the full 
five-dimensional space of parameters is not practicable and, 
as will be clear from the analysis below, is not warranted by 
the data. We neglect the intrinsic size of the pulsar emission 
region {Tt — >■ oo), and for PKS1257-326 we consider only the 
case Tb = 10^^ K (in the observer frame) as representative of 
the brightness temperature of a compact synchrotron source. 
We are also able to fix the scattering strength, q, for each 
set Qi, Qo) from the basic requirement on the image size 
or variability amplitude, which proved to be relatively stable 
from one realisation to another. Of the remaining three pa- 
rameters, the inner and outer scales have a restricted scope 
in that the theoretical appeal of a power-law is its scale-free 
nature. Our initial approach is therefore to set the inner 
scale to below our numerical resolution (10* cm or 10^ cm 
for PSR or IDV simulations, respectively), the outer scale 
to above the extent of the simulated domain (10^^ cm or 
10^^ cm), and thus proceed to constrain p. Subsequently we 
allow some flexibility in either Qi or Qo, as appropriate, to 
see how our constraints relax. 



3 PSR B0834+06 CONSTRAINTS 

The pulsar B0834-f06 offers one of the best examples of 
the parabolic arc phenomenon ( [Stinebring et aL]|2001[ |Hill| 
et al.|2005 l. The phenomenon itself is understood as arising 
from a quadratic relationship between the geometric delay 
and the scattering angle, for a preferred direction on the 
sky (Cordes et al. 20061. That preferred direction may be 
determined by the effective transverse velocity vector, but 
in cases such as B0834-I-06 where the arcs are very sharply 
defined one infers that the preferred direction arises from 



a strong anisotropy in the scattering (Walker et al. 2004 



Cordes et al. 20061. This inference has been confirmed by 



construction of the electric field image in the delay-Doppler 
plane (Walker et al 2008), and by intcrfcromctric imaging 
using very long terrestrial baselines (Briskcn ct al. 2010). 
By virtue of acquiring data simultaneously from several 
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Figure 2. Simulated scattered images of PSR B0834+06 for 
scale-free spectra with spectral indices (top to bottom): [5 = 
2, 2.5, 3, 3.5, 4, 4.5. Each image has a constant intensity of 10~^ 
added, to mimic the instrumental noise-floor. Tthe images are 
offset from each other by four decades along the vertical axis. 



sites, 



Brisken et al. (20101 were able to break some of the 



degeneracies which are present in single-dish data and thus 
determine, for example, the distance of the scattering mate- 
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rial along the line-of-sight {Dscr 
measured screen distance. 



3.1 Data and simulations 



415 pc); we adopt their 



pa 



To constrain our model of the scattering medium we use 
the electric-field image of Walker et al. (20081. For the 



present purposes, where we're modelling the scattering as 
one-dimensional, the appropriate mapping of the delay- 
Doppler image onto one angular coordinate is trivial: one 
simply tallies the squared magnitudes of the scattered waves 
at a given Doppler-shift. The resulting image is shown in fig- 
ure 1. With different physical approximations one can arrive 
at a different angular structure from the same delay-Doppler 



image; thus Gao et al.^ ( ,2010j and Rickett (20111 have ob- 
tained a two-dimensional rendering under the assumption 
that dispersive delays are negligible. Even so, the major part 
of their resulting image is approximately one-dimensional. 

The data shown in figure 1 form a single example of a 
scattered image; any hypothesised spectrum for the scatter- 
ing medium must plausibly be able to reproduce this image. 
It is not necessary to seek a direct match between specific 
features in the data and in the simulations. Rather we seek 
a match between the general characteristics of the observed 
and simulated images. Specific aspects of the simulated im- 
ages which we have paid attention to are: shape of the max- 
imum (plateau/peak/bell-like etc.); shape and slope of the 
profile "wings"; and the roughness of the profile (continu- 
ous/smooth/spiky etc.). We used the width of the image or 
its main component to fix a for each set (/3, Qi, Qo). 



3.2 Constraints from image profiles 

We first compared the data to simulated images for pure 
power-law inhomogeneity spectra (where the inner scale is 
unresolved and the outer scale is larger than the size of our 
simulated medium). Figure [2] shows a series of images repre- 
sentative of simulations for various values of the spectral in- 
dex, /3, with Qq^ = 10^^ cm and a negligible inner scale. The 
shallowest slopes, /3 < 2.5 produce image profiles that are 
too "triangular" , with the intensity trending steeply down- 
ward as one moves away from the origin. By contrast the 
observed profile (figure 1) exhibits multiple scales with, for 
example, a relatively flat region of the proflle terminating 
abruptly at 9 ~ —10 mas. 

Steep slopes, /3 > 3 are not compatible with the data ei- 
ther - the images they generate are too spiky or altogether 
fragmented in comparison with a fairly solid core to the 
observed image. Slopes 2.5 < /3 < 3 appear generally com- 
patible with the data. Note that we are simulating images 



in the 'snapshot' regime of the Narayan & Goodman ( 1989 1 
classification, where the telescope integration time is shorter 
than the diffractive time scale; this is appropriate for the 
data which we are comparing to. Previous work has shown 
that the data do not resemble the "ensemble-average" image 
expected for a Kolmogorov spectrum ( jBrisken et al.|[2bl0[ ). 

When the scale-free assumption is releixed, both shal- 
lower and steeper slopes can be made compatible with the 
data by restricting the range of the power law with a larger 
inner or a smaller outer scale, respectively. In particular, 
13 — 2 requires that > lO" — 10^° cm. For steep spectra 
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Figure 3. Electric field magnitudes of the scattered radiation 
in delay vs. Doppler shift coordinates. The upper panel is the 
Walker et al. (2008) construction from B0834+06 data, and the 
lower panel is a simulation with /3 = 3, Qq^ = 10^^ cm. For 
the B0834-I-06 data the locus of the power distribution conforms 
closely to a single parabolic arc, defined by the scattering geom- 
etry, whereas the simulation shows wiggles relative to the geo- 
metric arc. The wiggles are due to wave-speed (dispersive) de- 
lay contributions, and their absence in the data constrains the 
outer scale for steep fluctuation spectra. The intensity scale at 
the top is energy units relative to the maximum of a panel, 
10 log \u{t, fo)\^/ max \u{t, }u)\^- 



we find that /3 — 3.5 calls for Qq ^ < 10^^ — 10^^ cm, whereas 
for /3 — 4 a.n outer scale as low as Qq^ ~ 10^" cm needs to 
be chosen for simulations to resemble the data. 

We caution that the physical interpretation of our 
model (eq.[l| becomes more obscure as the inner scale, , 
increases and/or the outer scale, Qq^ , decreases. For a mean- 
ingful power-law model of the density fiuctuations we cer- 
tainly require ^ Q^^ - corresponding to the inertial 
range in the Kolmogorov theory, for example - and this is 
satisfied for all our models. But we note that for j3 > 3 
models to be consistent with figure 1 the outer scale must 
be comparable with, or even below the Fresnel-scale, which 
is If — 10^^ cm. Thus, although /3 > 3 fiuctuation spectra 
can be made consistent with figure 1, consistency can only 
be achieved by fundamentally altering the character of the 
spectrum to the point where it is no-longer steep in respect 
of its refractive {Q~^ > If) properties. The merit of /3 > 3 
models is therefore questionable. 
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3.3 Additional constraints from dispersive delays 

For models that produce acceptable image profiles, we also 
performed full dynamic-spectral simulations of the electric 
field. In some cases this revealed substantial deviations away 
from the notional parabolic locus of the power distribution in 
the delay-Doppler plane, as illustrated in the bottom panel 
of figure [3] We attribute these "wiggles" to dispersive (wave- 
speed) delay contributions. The absence of any clear wiggles 
in the reconstruction of Walker et al. ( 2008 1 (shown in the 
top panel of fig. [3| allows us to constrain the outer scale to 
be Qq ^ < 10^^ cm for 13 — 3. For steeper fluctuation spectra 
the outer scale is already limited to very small values by 
matching to the observed image profile. For Kolmogorov, or 
shallower slopes, there is relatively little fluctuation power 
at low spatial frequencies and for these spectra our delay- 
Doppler analysis added no new information on the outer 
scale. 

At this point it might be helpful to remind the reader that 
the specific inferences made in this section are restricted in 
scope to strictly one-dimensional models considered in this 
paper. 



4 IDV LIGHT-CURVE CONSTRAINTS 

It is possible that the plasma inhomogenities responsible for 
pulsar parabolic arcs are of the same basic nature as those 
responsible for Intra-Day Variability, as both phenomena 
appear to involve highly concentrated regions of very strong, 
highly anisotropic scattering material. In §5 we address that 
possibility. First, however, we consider what constraints can 
be placed on our one-dimensional power-law model (eq. [T]) 
from the IDV data alone. 



4.1 Data and simulations 



We utilised data for PKS1257-326 CBi gnall et "aL][2003 |. 
This source is one of three bright, compact radio quasars 
which are known to exhibit large-amplitude radio vari- 
ations on time-scales ~ hours. (The other two sources 
are PKS0405-385 [Kedziora-Chudczer et al 1997] and 
J1819-)-3845 [Dennett-Thorpe and de Bruyn 2000].) The 
data we employ are simultaneous dual frequency (4.8 GIfz 
and 8.64 GHz) light curves obtained with the Australia Tele- 



scope Compact Array (ATCA) as described in Bignall et al. 
(20031, supplemented by similar data taken in 2003 as part 



of the same observational program. 

We model the light curves assuming the Fresnel length 
of Ip — 3.3 X lO'' cm (in the X-band, as rescaled from the 



kinematical measurement of Walker et al. 20091, implying 



the distance to the scatterer of D ~ 6.4 pc. For simplic- 
ity, as the properties of the source are unknown in detail, 
we adopt equal angular sizes for the source in both C- and 
X-bands. (Assuming equal brightness temperatures instead 
does not affect the qualitative results derived below.) We 
use the observed amplitude of flux variation to fix a for 
each combination of {l3,Qi,Qo). 

Following the approach adopted in the previous section, 
we did not try to match simulations to the data exactly, fo- 
cusing instead on qualitative aspects of light curves. These 
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Figure 4. Correlation coefficient of PKS 1257-326 C- and X-band 
light curves and normalised inter-band lag, f for three observa- 
tional seasons. Positive lag means variations in C-band leading. 
For markers with error bars, the lag and correlation shown are 
averages for light curves from two or three consecutive dates. 



are best summarised by the adjectives "smooth" and "well- 
correlated", which our simulated light-curves must repro- 
duce. For source brightness temperatures Tt < 10^^ K, as 
is expected for synchrotron radiation, the size of the source 
has a strong smoothing effect on the scintillation, and as a 
result the smoothness of the light-curves tells us little about 
the phase fluctuation spectrum (for Tt in excess of 10^'' K 
the data strongly disfavour shallow spectra with inner scale 
below the Fresnel length). The degree and lag of inter-band 
correlation, on the contrary, prove useful to constrain the 
phase fluctuations at larger scales as the data exhibit very 
little prismatic shift. 



4.2 Constraints on prismatic shift 

The inter-band correlation statistics are summarised in Fig- 
ure [4] which shows the light curve pairs' correlation coef- 
ficient and inter-band lag normalisecQ by the overall vari- 
ability timescale: r = maxi- CCFx/c(''")/('Tix?ric) and f = 
argmaxr CCFx/c(^)/\/HWHMxHWHMc, where and 
HWHM are the height of the light curve ACF and half its 
width at the m^/2 level. As one can see, the light curves 
are strongly correlated in the observed regime of moder- 
ately weak scattering, with the coefficient rarely dropping 
below 0.7 and the lag being substantially below the tempo- 
ral extent of a typical "scintle" - i.e., the characteristic time 
scale of the observed quasi-sinusoidal variations of the flux 
density, quantified by the HWHM. 

The smallness of the observed lag strongly constrains 
the large-scale gradients, Vke, of the electron column den- 
sity, which are expected to induce prismatic shifts due to 
dispersion, i.e. dependence of the defiection angle, 9, on 
the wavelength. At the observer, the difference in arrival 
angles translates into a lag of tx/c ~ Ax/c^-Dcff/t'cff = 
(Aq — A^)feVKeDcff/27ri;ofi. Since all gradient modes of scale 
greater than a few light-curve lengths, T, contribute (inco- 
herently, hence quadratures) to the prismatic shift, the typ- 



^ Normalisation is introduced to take into account considerable 
annual variation in the effective velocity of the scattering screen. 
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ical lag, for a normal scale hierarchy Q^^ VcffT Qq 
is expected to be 



TX/C 



e, 



(tFlT)- 



/3 < 3 
/3 > 3 



(2) 



with the coefficient ~ and tp = Ip/vcS- 

For PKS1257-326 at C- and X-band, the scattering is 



moderately weak, implying a 



1/2 



0.1 and tp ~ du- 



ration of a scintle. Therefore for the prismatic shift to be 
small in Fresnel units we conclude from equation 2 that ei- 
ther (i) the slope of the fluctuation spectrum is shallow, 
or (ii) the slope is steep, but the outer scale is very small, 
Qo^ i$ ^-F ~ a fewx 10® cm. These conclusions are supported 
by our numerical simulations and are valid for various values 
of the source size around our assumed Tt ~ 10^^ K. 

We note that the observed lags have previously been 
interpreted in terms of an offset in the apparent centroids 
of the source in the C- and X-bands (Bignall & Hodgson 



2012[ [Bignall et al.|2003l Macquart et al. 2012).In this case, 
the measured values of tx/c should themselves be consid- 
ered as upper limits on the prismatic shifts introduced by 
the scattering screen, and the constraints we have given are 
therefore conservative. 



5 SUMMARY OF CONSTRAINTS 

We have constrained the form of the fluctuation spectrum 
(eq.[l} independently for PSR B0834+06 and IDV PKS1257- 
326. However, there is some commonality between the two 
sets of constraints: neither spectrum is permitted to have 
much refractive fluctuation — i.e. not much power in the re- 
gion IpQ < 1. In the scale-free case (Qilp — >■ oo, Qolp — > 0) 
that translates to /3 < 3 for both scattering media. This 
constraint on P relaxes if we admit a flnite outer scale, but 
for both pulsar and IDV source the outer scale is forced to 
lie close to the Fresnel-scale in order to expand the permit- 
ted range of /3 significantly on the steep side. We have al- 
ready noted the questionable merit of such a description for 
PSR B0834-I-06, and with the same circumstance required 
for IDV PKS1257-326 we reject the model. The reason is 
simply that the Fresnel-scale for our IDV data is two orders 
of magnitude smaller than for our pulsar data: steep spectra 
can only be made to fit the data by contrivance. 

Implicit in the above logic is the expectation that we 
can explain both PKS1257-326 and PSR B0834-f06 with a 
single fluctuation spectrum; we now concentrate on that is- 
sue. The known points of similarity are as follows. First, 
both require very strongly scattering media. Second, both 
require high levels of anisotropy in the scattering media — 
presumably enforced by a strong and ordered magnetic fleld. 
Third, there appears to be a high degree of spatial localisa- 
tion for both phenomena, in that there are no known cases 
in which a given scatterer heis been seen to affect more than 
a single background radio source. Fourth, as we have shown 
in sections 3 and 4, under the assumption of power-law fluc- 
tuations, both spectra are shallow, with /3 < 3. And finally, 
the recent measurement (Brisken et al 2010) of the parallel 
(to the scattering angle) velocity component of the scatter- 
ing medium towards PSR B0834-f 06 is similar to the results 



previously obtain ed for PKS1257-326 (|BignaU et al.||2003 l 
and J1819+3845 ( [Dennett-Thorpe fc de Bruyn|2003j . 

Our models allow us to reinforce the first of 
these points. The scattering strength inferred for 
PSR B0834-I-06, when rescaled to the X-band and 
the distance of the PKS1257-326 scattering medium, 

becomes ao834(Al257/Ao834)^(^F,1257/^F,0834)''~'^ — 

5.6 X 10^(321/8640)^(3.3 x 10V2.3 x lO")""^ ~ 0.016 
(numerical value for the scale-free /3 = 3 case). This is 
comparable to the actual Q1257 — 0.03 we required to model 
the observed scintillations of PKS 1257-326. The two values 
can be considered consistent with each other within the 
accuracy that the present analysis and data allow. In other 
words: both scattering phenomena can be attributed to the 
same type of physical structure, just placed at different 
distances from the observer. 

Given these points of similarity between the two scat- 
tering media, we consider it established that the scatter- 
ers are fundamentally similar. Thus our preferred spectral 
index range for power-law inhomogeneities is 2.5 < /3 < 
3; this includes the one-dimensional Kolmogorov spectrum 
{P = 8/3). At the uppermost end of this range the observed 
lack of wiggles in the B0834-I-06 parabolic arc constrains 
the outer scale to be Qq^ < 10^"^ cm. But for Kolmogorov, 
or shallower spectra these wiggles are too small to notice 
even when the outer scale is as large as permitted by inter- 
mittency of the IDV phenomenon (i.e. Qq^ < 10^* cm). 

Slopes shallower than /3 = 2.5 can be allowed by postu- 
lating a large inner scale, and our modelling does not rule 
out such possibility. However, similar comments apply in 
this case as to the interpretation of a steep spectrum with a 
very small outer scale. If a large inner scale is needed then 
we should consider the possibility that our model is inap- 
propriate, even though it can yield results similar to the 
observations. For example, it has been proposed by [Brisken] 
et al.[([2010[) that for B0834-I-06 the statistics of the electron- 



density fluctuations are non-uniform. If so then a great deal 
of freedom is permitted in the modelling, and even if the fluc- 
tuations are "power-law" none of that freedom is present in 
the model we have utilised here. 



6 NATURE OF THE SCATTERING MEDIA 

Unlike the strength parameter, a, the variance of the elec- 
tron column density is not fixed by the data, because the 
integral of ([T]) is dominated by its lower limit and the outer 
scale is poorly constrained for our preferred, shallow slopes. 
Luckily, this is not the case for the variance of the electron 
volume density. Assuming the overall screen thickness Z is 
of the same order as the outer scale, the electron volume 
density variance, for a generic orientation of the screen with 
respect to the line of sight, is 



io 2 



,3-/3 



7t[1 + (Qo/Q,)] aQi 
2r(^)sinf(/3-l) X^ril^r'' 



(3) 



which is very weakly dependent on Qo for /3 close to 3. In 
particular, for the scale-free P — 3 model the implied elec- 
tron density r.m.s. is Aue ~ 30cm~^, independent of the 
outer scale. And for the one-dimensional Kolmogorov spec- 
trum this reduces to Arte ~ 10 cm~^ for the largest plausible 
outer scale (10^* cm). This value is consistent with Rickett's 
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(2011) estimate for B0834+06, based on the same data we 
have used. 

Bearing in mind that we are dealing with ionised gas 
(temperature ~ 10* K), the foregoing estimates imply that 
the scattering medium manifests fluctuations in gas pressure 
which are an order of magnitude greater than typical pres- 
sures in the diffuse ISM: Pgas 3> Pism- This is puzzling, as 
one expects pressure equilibration on a sound-crossing time. 
This is a familiar problem in the context of Extreme Scat- 
tering Events (Romani, Blandford & Cordes 1987), where 
the implied gas pressures are even larger than our estimate 
following equation (3). But for the scattering media we have 
modelled, assuming that the anisotropy is due to a magnetic 
field, the gas pressure problem is only the tip of the iceberg. 

The high anisotropy in the electron-density fluctuations 
indicates that the gas pressure fluctuations are small in com- 
parison with the magnetic stress in the medium: Pb ^ Pgas- 
And the magnetic held must have its origins in a current sys- 
tem which is held in place by material stresses. The speciflc 
nature of the material - call it "X" - is unknown; but in or- 
der for the material stresses to at least balance the magnetic 
stresses we must have Px ^ Pb- Thus we have 



Px>Pb^ Pnas > Pi. 



(4) 



This ordering indicates that the diffuse ISM has little in- 
fluence upon the scattering media we are modelling, and 
is practically irrelevant to the dynamics of the material, X, 
which is ultimately responsible for the observed phenomena. 

The conclusions we have just drawn rest on the assump- 
tion that Z ~ Qo^ - If, instead, Z is much larger than the 
outer scale then smaller pressure fluctuations correspond to 
a given column-density variance. However, it is only possi- 
ble to mitigate the stress hierarchy of equation 4, so that 
Px ~ Pism, ii Z > 10^ Q^'^ — i.e. an enormous geomet- 
ric factor must be invoked. A further, more serious prob- 
lem faces us if we try to overthrow the stress hierarchy of 
equation 4. If Pb ~ Pism, so that the magnetic stresses 
are no higher than typical interstellar stresses, then we do 
not expect the magnetic held to have a constant orientation 
through the scattering medium. But the fact that the ob- 
served anisotropy is A ^ 1 tells us that the orientation is 
constant to within a projected angle of < 1/A radian. For 
these reasons we consider "ambient pressure" models of the 
IDV and parabolic arc phenomena (e.g. Pen fc King||2012 l 
to be untenable. 



7 SUMMARY AND CONCLUSIONS 

We have numerically simulated radio-wave scattering by 
a totally anisotropic, power-law spectrum of electron- 
density fluctuations. Comparing our results to data on PSR 
B0834-I-06 and IDV PKS1257-326 we constrain the spectral 
index of the power-law to lie in the range 2.5</3<3 — a 
range which includes the one-dimensional Kolmogorov spec- 
trum. The permitted range of /3 can be expanded by intro- 
ducing inner and/or outer scales, but signiflcant expansion 
requires such a large inner scale, and/or small outer scale, 
that one or both scales are playing a major role, thus de- 
stroying the essence of the power-law. 

We flnd that the amplitude of the density fluctuations 
is comparable in the two distinct scattering regions which 



we model. This, together with the small spatial extent of 
the scattering media and their high levels of anisotropy, en- 
courages the view that some pulsar parabolic arcs and the 
IDV of flat-spectrum quasars are due to the same types of 
scattering concentrations. The nature of these scatterers is 
obscure, but a simple-minded consideration of the associ- 
ated physical stresses shows that there is a strong hierarchy, 
with the diffuse ISM at the lower end. We conclude that the 
diffuse ISM is largely irrelevant to these remarkable scintil- 
lation phenomena. 
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